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We investigate the propagation of partially coherent beams in spatially nonlocal nonlinear media 
with a logarithmic type of nonlinearity. We derive analytical formulas for the evolution of the beam 
parameters and conditions for the formation of nonlocal incoherent solitons. 
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Incoherent solitons are partially coherent optical 
beams propagating without changing their shape in ma- 
terials with a slow nonlinear response 0,0- Such nonlin- 
ear materials respond to the time averaged light intensity, 
which means that even though the phase of the partially 
coherent light beam fluctuates randomly, it still leads to a 
smooth light-induced refractive index change and subse- 
quent trapping of the beam dSSHHEIa. 

Incoherent 

solitons have been experimentally observed almost solely 
in photorefractive media, which exhibit a strong and suf- 
ficiently slow saturable Kerr- like nonlinearity 0, Q . 

In inertial bulk Kerr media (2+l)-dimensional par- 
tially coherent beams are still unstable and will either 
diffract or self-focus unboundedly, dependin g on whether 
the power is above a certain critical value |T(i ITU Il2| . 
Incoherence counteracts self-focusing and leads to an in- 
crease of the critical power, but it cannot remove the col- 
lapse instability [id 1 1 1 1 ] . Consequently all experimental 
observations of incoherent solitons have been in saturable 
materials, such as photorefractive materials and liquid 
crystals 

Recently Peccanti et al. reported the first obser- 
vation of incoherent solitons and their interaction in 
nematic liquid crystals |13| . Liquid crystals have a 
strong non-instantaneous saturable Kerr-like nonlinear- 
ity associated with light induced molecular reorientation 
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Their nonlinear res pon se time ranges 
from tens to hundreds of milliseconds [Tj| , which is suffi- 
ciently slow to allow for the formation of incoherent soli- 
tons. Apart from being inertial the nonlinear response 
of liquid crystals is also inherently spatially nonlocal, 
because the molecular reorientation induced by a light 
beam in a particular place will affect the orientation of 
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molecules far beyond this point, due to the long-range 
character of their interaction [lj, [lg] • 

Spatially nonlocal nonlinearities are in fact an inher- 
ent property of many physical systems, including matter 
waves, nonlinear optics, and field theory. As such the 
nonlocal aspect of nonlinear interaction has attracted sig- 
nificant interest recently [H^ . I20I ] . Nonlocality can have 
profound consequences on the properties of optical beams 
and soliton formation, e.g., leading to collapse arrest of 
finite-size beams 0, attraction and formation of bound 
states of dark solitons [2^], and modulational instabil- 
ity (MI) in defocusing materials [2^, [24|. A nonlocal 
nonlinearity may even describe parametric wave mixing 
and solitons j2^. In the particular case of liquid crys- 
tals certain aspects of their nonlocal nonlinear response 
have already been discussed in the literature, including 
theoretical and experimental studies of the formation of 
multiple solitons [2H. |27| , accessible solitons [2£| as well 
as soliton interaction and MI |2|j . 

Motivated by the documented unique dual inertial and 
nonlocal character of the nonlinearity of liquid crystals 
and the recent experimental observation of incoherent 
solitons in this material, we analyze here the effect of non- 
locality on the propagation of partially coherent beams 
and formation of incoherent solitons. Because such a 
problem is in general analytically intractable we consider 
here the special case of a nonlinearity with a logarith- 
mic dependence of the refractive index change on light 
intensity. This model has previouslybeen used to study 
incoherent beams in local media 0,H H3, El [13 as weu 
as fully coherent beams in nonlocal media |33| . Here 
we extend it to describe simultaneously incoherence and 
nonlocality. While being rather specific, the logarithmic 
model is unique in the extent that it enables a fully ana- 
lytical treatment of the nonlinear beam propagation. It 
has been successfully used to predict a variety of effects 
associated with beam propagation and soliton formation 
ii. 
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Propagation of a two-dimensional quasi- 
monochromatic partially coherent beam with the 
slowly varying amplitude ipif, z) is governed by the 
following nonlinear Schrodinger equation: 
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(1) 



where r=(x,y), I=\tp\ 2 denotes the light intensity, Sn(I) 
represents the light induced refractive index change, and 
Vl = d 2 /8x 2 + d 2 /dy 2 . 

In the following we use the mutual coherence function 
rfn , ry, z) to represent the coherence properties of the 
beam 35] . This function describes the spatial correlation 
between the field at points f\ and r*2, and is defined as 



(2) 



where brackets denote temporal or ensemble averaging. 
In terms of the mutual coherence function the (time av- 
eraged) intensity is given by 



I{f, z) = r(r, r, z). 



(3) 



We take the refractive index change Sn to be the following 
nonlocal function of light intensity 



5n(f, I) = ri2 In 



(4) 



where J dr = J_ J_ oa dxdy and R(r)=R(r) is the so- 
called nonlocal response function, which depends only on 
the length of its argument vector, r = |rj = y 1 x 1 + y 2 . 

The model (0J is phenomenological and reflects the fact 
that the refractive index change in a particular point in 
space is not only determined by the wave intensity in this 
point, but also on the intensity in a certain surround- 
ing region given by the width of the response function. 
The width of the nonlocal response function therefore 
determines the degree of nonlocality. In the special case 
when R(r) = <5(r) this model describes a local nonlinear 
medium 5n = ri2m/. A simplified model may also be 
derived in the weakly nonlocal limit, in which the width 
of the response function is finite, but small compared to 
the beam width [2(],[36j]. In the strongly nonlocal limit, 
in which the response function is much broader than the 
beam, the governing model becomes linear |20t l33| . 

In general, the particular form of the response function 
is determined by the specifics of the physical process re- 
sponsible for the nonhnearity of the optical medium. For 
instance, it can be shown that for the reorientational non- 
linearity of liquid crystals and general diffusion type non- 
linearities, the nonlocal response can be approximated as 



R(r) oc exp(— r/a), 



(5) 



with a representing the extent of the nonlocality. How- 
ever, the majority of nonlocality-mediated effects appear 
to be rather generic and do not depend strongly 

on the particular form of the nonlocal response function. 



Here, to enable us to make an analytical description, we 
therefore assume that the nonlocality is described by the 
normalized Gaussian response function 



1 ( x z + 
R(r) = — - exp - 

7T(T Z \ (7 
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While we consider here a nonlocally isotropic medium, 
the model can be easily extended to include anisotropy 
of the nonlocal properties, such as those caused by a di- 
rectional heat flow in thermal media j37| ■ 

The mutual coherence function satisfies the following 
propagation equation 0, ^| 
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Introducing the independent spatial variables p and q, 
P=~(ri+f 2 ), q = fi-f2, (8) 

and using the specific form of the nonhnearity trans- 
forms the propagation equation into 

dT 

^ + v p -- v,j-r + 

\jR(f 2 -mm 



(9) 



In what follows we will consider the so-called Gaussian- 
Schell model for the incident partially coherent beam |35| . 
We will also restrict our analysis to circular beams, i.e. 
beams having isotropic coherence properties. The ex- 
tension of the analysis to elliptical incoherent beams is 
straightforward. The mutual coherence function of the 
Gaussian-Schell beam is expressed as 



r(p, q, z = 0) = exp 



Pl 



(10) 



where the initial effective coherence radius #o is given by 
the relation 



l/6l = l/r 2 c + l/(4p 2 ). 



(11) 



Here po is the initial radius and r c is the initial coher- 
ence radius of the beam, respectively. Due to the loga- 
rithmic form of the nonhnearity, the beam will maintain 
the Gaussian statistics during propagation and thus the 
coherence function will keep the form given by Eq. i|10|) 
We therefore look for the solutions to Eq. © using the 
following Gaussian ansatz 



r(p, q, z) 
A(z) exp 



lit 
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(12) 



where A(z) and p(z) represent the amplitude and phase 
variations of the coherence function, and p(z) and 0(z) 
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are its diameter and generalized coherence radius, respec- 
tively. The input amplitude plays no role for the dy- 
namics and thus we use the initial conditions A(Q)=1, 
p(0)=po, and 8(0)=8o, Inserting the ansatz l(T2l into 
Eq. © leads to a set of ordinary differential equations 
for the parameters of the coherence function 
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(13) 
(14) 
(15) 
(16) 



From the first two equations i|13|) and i|14|) we obtain 
the relation 



0/p= Qo/Po, 



(17) 



which shows that during the evolution the beam con- 
serves its coherence, defined as the number of speckles 
within the beam diameter. Interestingly, this relation 
does not depend on whether the nonlinearity is local or 
not. Combining Eq.{nj, Eq.JTJJl and Eq.JT2J we obtain 
the following evolution equation for the beam radius 



fp 

dz 2 



A 



2n 2 



P 



(18) 



This equation describes the dynamics of the width of 
a partially coherent beam with Gaussian statistics in a 
nonlocal medium with a logarithmic nonlinear response. 
Taking a = we recover the expressions governing in- 
coherent beams in a local medium |3jJ. On the other 
hand, for r c — > oo, Eq. Ijl8|l describes coherent beams in 
a nonlocal medium. It is clear that the main role of the 
nonlocality is to effectively weaken the influence of the 
nonlinearity on the propagation of the beam. For very 
high degrees of nonlocality the nonlinear effects becomes 
negligible and the beam just diffracts. 

Incoherent nonlocal soliton solutions are obtained from 
Eq. H18|) by setting p(z) = po. This gives the relation 



n 2 



1 



Pi 



(19) 



As expected, the existence of bright incoherent solitons, 
or localized beams, is seen to require a focusing nonlocal 
nonlinearity with n 2 > 0. Importantly, the relation i|19|) 
further shows that trapping of a beam with a given initial 
width po and coherence radius r c in a nonlocal medium 
requires a stronger nonlinearity (higher n 2 ) than in the 
case of a purely local nonlinear response. This is again 
due to the fact that the nonlocality effectively leads to a 
decrease of the strength of focusing and subsequently to 
a weaker localisation of the beam. From Eq. i|19|) , we find 



the expression for the radius of the incoherent nonlocal 
soliton 



2 Po 2 
Po = — 



4a 2 
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PO 2 



where 



Po 2 = Po{& = 0) = 



2n 2 - A/r 2 c 



(20) 



(21) 



denotes the radius of the incoherent soliton in a local 
medium. 
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FIG. 1: (a) Soliton width (p) as a function of the degree of 
nonlocality (a) and coherence radius (b). Here n2 = 1.8 

In Fig. H^a) we plot the radius of incoherent nonlocal 
solitons versus the degree of nonlocality a for three values 
of the coherence radius r c . These plots clearly show that 
the soliton width increases with nonlocality and becomes 
proportional a in the highly nonlocal regime. One can 
show that for a/r c 3> 1 



Po = 2p — . 

r r 



(22) 



The degree of nonlocality is an inherent property of the 
medium and as such difficult to control. However, one 
can easily vary the coherence properties of an incident 
beam. Therefore we plot in Fig. ^b) the soliton width 
po versus the coherence radius r c for three values of the 
nonlocality a. Again, the increase of the soliton width 
with incoherence and nonlocality is evident. Addition- 
ally, Fig. njb) shows that for a given strength of nonlin- 
earity there is minimum coherence radius, below which 
the incoherent soliton cannot exist. The existence of such 
a threshold is well-known from earlier works on incoher- 
ent solitons in local media 0, 0, |^|. Interestingly, it 
turns out that for a nonlocal nonlinearity, this threshold 
is also determined solely by the strength of the nonlin- 
earity and the degree of coherence through the relation 



n 2 - 2/r 2 c > 0. 



(23) 



To understand why this condition does not depend on 
the beam size or the degree of nonlocality we recall that 
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the soliton is formed when the nonlinearity fully com- 
pensates the transverse spreading of the beam, which is 
caused by incoherence and diffraction. The former is due 
to the diffusive character of the beam, while the latter is 
a direct consequence of the finite size of the beam. For 
a very broad beam the contribution to its expansion due 
to diffraction and nonlocality becomes negligible com- 
pared to that due to incoherence. Hence, to achieve beam 
trapping the strength of the nonlinear response of the 
medium should be high enough to compensate at least 
the incoherence-induced beam spreading. This, in fact, 
determines the threshold for soliton generation. As soon 
as this condition is satisfied a bright incoherent soliton 
can exist, whose width will be determined by the inter- 
play of nonlinearity, coherence, diffraction, and nonlocal- 
ity. Obviously, the resulting beam diameter will always 
be larger than that of the local soliton of the same degree 
of coherence. 

More physical insight into the competition between 
diffraction, nonlinearity, nonlocality, and incoherence can 
be obtained by an effective particle analogy. Assuming 
that (dp/dz)(z—0)=0, and integrating Eq. l |T8|l once, one 
obtains a classical-mechanical equation 



{dp/dzf + P(p) = 0, 



(24) 



describing an effective particle p with kinetic energy 
(dp/dz) 2 moving in the potential P(p), given by 



P(P) 



ff 2 



- 1 



2?i9 In 



Pi 



(25) 



The asymmetric potential is depicted in Fig.|3for a beam 
with effective coherence radius 9q— 1 moving in a focusing 
logarithmic nonlinear medium with 72.2=1.8. As nonlocal- 
ity increases the width of the potential also increases and 
its minimum becomes less pronounced. 
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FIG. 2: (a) Potential P(p) in a focusing logarithmic nonlinear 
medium (ri2=l-8) for different degrees of nonlocality a. The 
initial beam parameters are 9q — 1 and po = 1. 

Stationary soliton solutions with constant beam width 
correspond to the effective particle being located at the 



bottom of the potential well. For initial beam parame- 
ters different from the soliton solution the beam width 
(and coherence radius and peak intensity) will undergo 
periodic oscillations, corresponding to the effective parti- 
cle oscillating in the bottom of the potential well. As the 
profile of the potential well changes with higher a one can 
expect that the amplitude and period of beam oscillations 
will drastically increase for high degree of nonlocality. In 
the limit of the small amplitude oscillations, their form 
can be obtained analytically in the form of Jacobi ellip- 
tical function, in the way analogous to that discussed in 

m 

This is indeed the case as illustrated in Fig. |3| where 
we show the nonstationary propagation of the partially 
coherent beam in nonlinear medium. The initial parame- 
ters are chosen such that po = 1.0, r c = 1.15, and n-i = 1. 
In the graph we plot the beam radius as functions of the 
propagation distance in case of local (dotted line) and 
nonlocal (solid line) nonlinearity. The increased ampli- 
tude and period of oscillations induced by the nonlocality 
is evident. 
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FIG. 3: Illustration of nonstationary propagation of a par- 
tially coherent beam in a logarithmic nonlinear nonlocal 
medium. The nonlocality parameter varies as a = (dot- 
ted line) and a = 10 (solid line). 

In conclusion, we analyzed the propagation of incoher- 
ent optical beams with Gaussian statistics in a nonlinear 
nonlocal medium. Considering a special case of a log- 
arithmic type of nonlinearity and Gaussian nonlocal re- 
sponse we obtained an analytical formulas governing all 
relevant parameters of the beam. We showed that non- 
locality results in increased diameter of the incoherent 
solitons. When the initial parameters of the beam differ 
from the exact nonlocal soliton solution, beam experi- 
ences periodic expansions and contractions whose period 
increases with nonlocality. 
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